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ABSTRACT 

Context. Core-collapse supernovae (CC SNe), especially those of type Il-plateau (II-P), are thought to be important contributors to 
cosmic dust production. The most obvious indicator of newly-formed and/or pre-existing dust is the time-dependent mid-infrared 
(MIR) excess coming from the environment of SNe. In the past few years several CC SNe were monitored by the Spitzer Space 
Telescope in the nebular phase, hundreds of days after explosion. On the other hand, only a few of these objects have been analyzed 
and published to date. 

Aims. Our goal was to collect publicly available, previously unpublished measurements on type II-P (or peculiar IIP) SNe from the 
Spitzer database. The most important aspect was to find SNe observed with the Infrared Array Camera (IRAC) on at least two epochs. 
The temporal changes of the observed fluxes may be indicative of the underlying supernova, while spectral energy distribution (SED) 
fitting to the fluxes in different IRAC channels may reveal the physical parameters of the mid-IR radiation, which is presumably 
caused by warm dust. 

Methods. The IRS spectra were extracted and calibrated with SPICE, while photometric SEDs were assembled using IRAF and 
MOPEX. Calculated SEDs from observed fluxes were fit with simple dust models to obtain basic information on the dust presumed 
as the source of MIR radiation. 

Results. We found twelve SNe that satisfied the criterion above, observed at late-time epochs (typically after +300 days). In three 
cases we could not identify any point source at the SN position on late-time IRAC images. We found two SNe, 2005ad and 2005af, 
which likely have newly formed dust in their environment, while in the other seven cases the observed MIR flux may originate from 
pre-existing circumstellar or interstellar dust. Our results support the previous observational conclusions that warm new dust in the 
environment of SNe contributes only marginally to the cosmic dust content. 

Key words, supernovae: general ~ supernovae: individual - dust, extinction - infrared: stars 



1. Introduction 

According to theoretical and observational studies during the 
past four decades, dust grains may be relevant components of the 
circumstellar matter (CSM) in the environment of core-collapse 
supernovae (CC SNe), which are the last stages of the life cy- 
cles of massive (M > 8 Mg) stars. However, while the possibil- 
ity of dust fo rmation around CC SNe was first raised more than 
40 ye ars ago dCernuschi et al.L[l967l : lHovle & W ickrama singhel 
Il970h . there are still many open questions in this topic. In the 
local Universe, intense mass-loss processes of asymptotic gi- 
ant branch (AGB) stars are thought to be the primary sources 
of dust. At the same time, the large observable amount of dust in 
high-redshift galaxies is hard to explain with AGB stars, espe- 
cially because their formational timescale is longer than the age 
of these distant ga laxies (there are some deta i led reviews abo ut 
this topic, see e. g. lDwek & CherchneffllloTltlGall et an.l201 11) . 

This is the point where CC SNe become relevant as ma- 
jor sources of dust in the early Universe. The lifetimes of their 
progenitors are usually shorter than ~1 Gyr (there ar e dust- 
rich hi gh-z galax ies with this age, see e. g. Maioli no et al.L 
l2004tibwek et al.L l2007h. The classical mod e ls of dust produc- 
tion rates of CC SNe dKoz asa et all '1989"; 'Todini & FerraraL 
|2001t iNozawaeLali |2Q03t [Morgan & Edmunds, 2003) all sup- 
port this theory, predicting very similar masses (0.1-1 M©) of 
newly formed dust. In more recent studies Bianchi & Schneiden 
(2007); Nozawa et alj ( l2007h : iKozasa et al] ( f2009h : ISilvia et alj 
t2010.) also presented similar results. 



An important conclusion of the latter papers is that the 
amount and the basic properties of the dust formed in the ejecta 
depend strongly on the type of the proge nitor, especially the 
thickness of its outer H and/or He layers. As lKozasa et alj (l2009h 
showed, a thicker outer envelope results in a lower expansion 
velocity and higher gas density in the deeper layers that con- 
tain the condensable elements (C, O, Mg, Al, Si). Therefore, 
larger (>0.03 fim) grains are able to form. This in turn means 
a larger initial dust mass and also a higher survival rate of the 
grains. This rate is also increased by the delayed arrival of the 
reverse shock at the dust-forming region (because of the lower 
gas density in the shocked region). Based on these theoretical 
expectations, type II-P SNe are the best candidates for dust for- 
mat! on_am2ngCC_SNe. Calculations of dust production efficien- 
cies (iGall etalil20Tll) also show that the mass range of 8-20 Mq 
(which is thought to be the typical progenitor mass range of II-P 
SNe) is the most important for all SNe. 

However, grain formation in the expanding and cooling 
ejecta is not the only way to find dust around CC SNe. If there is 
a relatively dense CSM around the progenitor (as it is in type Iln 
SNe), it has probably been formed via mass-loss processes from 
the progenitor. This pre-existing matter also contains dust grains 
and may also permit grain condensat ion interacting with the for- 
ward shock (see more details e. g. in iFox etalll2010ll2011l) . 

Although the results of different models seem to be con- 
cordant, observations of CC SNe in the local Universe do not 
support the prediction of the intensive dust production in these 
objects. Optical and mid-infrared analysis of the famous SN 
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1987 A in the Lai-ge Magellanic Cloud (LMC) revealed -10 
Mq of newly formed dust in the ejecta a few hundred days af- 
ter explosion (see e. g. D anziger et al., 1989, 1991; Lucv et al., 
1 198911 199 It iRoche et all Il993t I Wooden et all 11993 1. which is 
orders of magnitude lower than the value predicted by theoreti- 
cal studies. Similar results were found in other CC SNe studied 
in the past two decades based primarily on the mid-lR data of 
the Spitzer Space Telescope, hereafter Spitzer. None of these ob- 
jects were unambigously found to have more than 0.01 dust 
in their environment (see Section|4]for more details). 

This controversy has been tried to be resolved in differ- 
ent ways. There is a widespread agreement that theoretical 
models are not perfect, and there are both micro- and macro- 
physical effects conn ected with dust formation that are not 
under stood (see e. g. ICherchneff & DweE 120101: iFallest etaTl 
1201 lb- Circumstel lar dust may form optically thick clumps 
dLucv et aLl [1989'), resulting in a dust mass higher by an or- 



der of magnitude in nume r ical calculations ( Sugerman et all 
l2006t lErcolano et all 120071: iFabbri et all l201lh . but it is still 
not a satisfying explanation. There is also a possibility that 
the missing dust mass is hidden in the form of cold (< 50 
K) grains, whose thermal radiation is mostly invisible for 
Spitzer'^, instruments. This cold dust should stay mainly in 
the SN ejecta sev eral years after explos i on, in the so - called 
transitional phase (Sugerm an et all 120121: iTanaka eFal I I20T2L 
[jemim et al., 2012a b), or in older supernova remnants (SNRs). 
Some of the nearby SNRs were observed at mid-lR, far-lR, 
and submillimeter (sub-mm) wavelengths, but the inferred dust 
masses vary on a large scale (from 10"^ to 1 Mp), sometime s 
even for the same object (see the review of lGallet all 1201 Ih . 
Recently, observations from the Herschel Space Observatory 
(hereafter Herschel) are expected to give reassuring answers. 
The published results suggest relatively large dust masses: 
0.08 M0 in Cas A (Barlow et al., 2010), 0.1-10 Mq in N49 
in the LMC (Otsuka et al.L |2010D, and 0.4-0.7 Mq in SN 
1987A (iMatsuu ra et al.,. .201 11) : although new measurements on 
SN 1987A made with the Atacama Pathfinder Experiment 
(APEX) at sub-mm wavelengths do not support the presence of 
such a large amount of dust (Lakicevic et al, 2012). 

There are also some other ideas, but they are not supported 
by any convincing observational evidence. The observed low 
dust production rates of CC SNe could be compensated with 
a top-heavy initial mass function (IMF), which might be a 
characteristic of some e arly galaxies ( e. g. Bromm e t al., 2002; 
iTumhnsonl |2006; Mic halowski et al 1 I2010r It is also possi- 
ble that pair-instability supernovae (PISNe), the assu med fate 
of s tars with progenitor masses of 140-260 M© (Hege r et all 
12003 ). could be e fficient dust pro ducers in the first galax- 
ies ( jNozawa et al.j 120031: iGall et all ,2011), but there are still 
no strong evidence for the existence of these objects. Some 
authors argue that depending on star-formation history, AGB 
stars may b e significant dust sources in galaxies younge r 
than 1 Gyr ( Vali ante et all 120091: iDwek & Cherchnefil l201lh . 
There is another intere sting idea about p ossible grain conden- 
sation in quasar winds (lElvis et all l2002h . Beyond these ideas, 
more calculations suggest significant grain growth in the ISM, 
which may be an acceptable explanation of both the low dust 
production rat es of SNe and the high dust content of dis- 
tant galaxies (iDraind. l2003[ l2009t iMichalowski eTall 1201 Ot 
'A 



iMattsson & Andersenfl2012l:iMattsson et all 1201 IT 



To get closer to the mystery of the origin of cosmic dust, 
it is necessary to have a large amount of high-quality observa- 
tional data. Up to now, there are only 20-25 SNe for which we 
have quantitative pieces of information about their dust content; 



see e. g. lFox et al ] (120 1 Ih . iGail et al.1 (l20 1 ih . ISzalai et al.1 (l20 1 ih . 
and Section |4] for details. Regarding the type 11-P SNe, there 
are less than ten objects with detailed analyses, and we know 
quite a few examples without detectable dust content. iFox et al.l 
(I2OI ll) published the results of a survey to detect MIR excess 
(the most obvious evidence of warm dust) in the environment 
of type lln SNe. They found late-time emission in ten objects 
from their sample of 69 SNe and concluded that the lln subclass 
represents the largest part of dust-containing SNe (it should be 
noted that the authors carried out this project during the warm 
mission of Spitzer, so they were able to use only the 3.6 and 4.5 
jim channels supplemented by optical and near-lR photometry 
and spectroscopy). 

Motivated partly by the work of lFox etaTI ( 1201 lb . our goal is 
to perform such analyses for type II-P SNe that are thought to be 
the best candidates for producing newly formed dust. We collect, 
reduce, and analyze unpublished measurements on CC SNe be- 
longing to the 11-P subclass, found in the public Spitzer database. 
In Section|2]we show the steps of data reduction and present the 
extracted photometric and spectroscopic data. We present ob- 
served spectral energy distributions (SEDs) of the studied SNe, 
and compare them with different dust models (Section O. We 
note that although the properties of the cold ISM are better con- 
strained by far-IR and/or sub-mm data, the low angular resolu- 
tion of the available instruments prevent measurements of ex- 
tragalactic SNe. The Spitzer data we applied in this paper probe 
only the warm (T ~ 200 - 1000 K) dust in the environment of ex- 
tragalactic SNe. Thus, our study focuses on the presence/absence 
of warm dust and its physical properties. We discuss the results 
of the SED fitting in Section [32] Finally, we summarize our re- 
sults and conlusions, and give a brief discussion on the role of 
type II-P SNe in cosmic dust production. 

2. Observations and data reduction 

We queried the public Spitzer database to find type II-P SNe 
observed with the Infrared Array Camera (IRAC) on at least 
two epochs. We found 12 SNe satisfying this criterion. We also 
collected observations obtained by Spitzer's Multiband Imaging 
Spectrometer (MIPS) and Infrared Spectrograph (IRS) instru- 
ments (not only the spectra but also the broad-band photomet- 
ric data recorded in peak-up imaging mode, PUI). In addition, 
we also checked the literature for any other observational data 
obtained around the same epoch as the Spitzer measurements. 

2.1. Supernova sample observed by Spitzer 

SN 2003 J was discovered in NGC4157on 2003 Janu ary 11.3 
UT (dates are in UT thereafter) bvlKushida et alJ (l2003h (~ 16.7 
mag, unfilteredV lAvani et al.l (120031) cl assified 2003J as a n ormal 
type-II SN. Using the infrared map of lSchlegel et alj (11998 ), the 
calculated value of galactic interstellar reddening is E{B—V)gai - 
0.021 mag. The estimated distance (D) to the host g alaxy (based 
on su rface brightness fluctuations) is 14.7 Mpc dTonrv et all 
I2001h . 

SN 2003hn was discovered in NGC 1448 on 2003 August 
25.7 by Evans et al. (2003) {V ~ 14.4 mag). Preliminary spec- 
troscopic analysis by Salvo et all (l2003l) showed that 2003hn 
was a type II-P supernova. Detailed analysis based on optical 
and near-inf rared photometry and optical sp ectroscopy was car- 
ried out by Krisciunas et al.J (120091) . while iJones et al.l (l2009l) 
and Olivares et al. (2010) also analyzed the data of SN 2003hn 
and of other type II-P SNe to obtain precise distances to the 
host galaxies. We adopt 18.1 Mpc as the distance of the SN 
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Table 1. Basic data of the supernovae 



Name 


RA 


DEC 


Galaxy 


Date of 


MJD - 


D (Mpc) 


E(B-V) 


References 




(J2000) 


(J2000) 




explosion 


2450000 




(mag) 




SN 2003J 


12:10:57.7 


+50:28:31.8 


NGC4157 


2003-01-11+ 


2651+ 


14.7 


0.021* 


1, 2 


SN 2003hn 


03:44:36.1 


-44:37:49.0 


NGC 1448 


2003-08-07 


2859 


18.1 


0.19 


3, 4,5 


SN 2003ie 


12:03:18.1 


+44:31:36.8 


NGC4051 


2003-08-16 


2868 


15.5 


0.013* 


1, 6, 7, 8 


SN 2004A 


16:43:01.9 


+36:50:12.5 


NGC 6207 


2004-01-06 


3011 


20.3 


0.06 


9 


SN 9nn'larl 






NGP Q41 




3407+ 








SN 2005af 


13:04:44.1 


-49:33:59.8 


NGC 4945 


2005-01-07 


3379 


3.9 


0.183* 


1, 11 


SN 2005CS 


13:29:52.8 


+47:10:36.1 


NGC 5194 


2005-06-28 


3550 


8.4 


0.04 


12, 13, 14, 15 


SN 2006bc 


07:21:16.5 


-68:59:57.3 


NGC 2397 


2006-03-24 


3819 


14.7 


0.205* 


1, 7, 16 


SN 2006bp 


11:53:55.7 


+52:21:09.4 


NGC 3953 


2006-04-07 


3833 


17.5 


0.400 


17 


SN 2006my 


12:43:40.7 


+ 16:23:14.1 


NGC 4651 


2006-07-26 


3943 


22.3 


0.027* 


1, 18, 19, 20 


SN 2OO60V 


12:21:55.3 


+04:29:16.7 


NGC 4303 


2006-08-16 


3964 


12.6 


0.022* 


1, 7, 19 


SN 2007OC 


22:56:41.8 


-36:46:22.3 


NGC 741 8A 


2007-11-03+ 


4408+ 


28.0 


0.023* 


1, 21 



References. (1) Schlegel et al. ( 1998); (2) Tonrv et al. (2001); (3) Jones et al. (2009); (4)'Krisciunas et al.' ('2009); (5) Olivares et al.' ("2010); (6) 
tierce &Tullv (1988);"(7) Smartt et al. (2009); (8) Harutvunvan et al. (2008); (9) Hendrv et al. (2006); (10) Li et al. (2011); (11) Kotak et all 
(I2OO6); (12) Pastorello et al. (2006); (13) Baron et al. (2007); (14) Pastorello et al. (2009); (15) Vinko et al. (2012); (16) Gall agher et aiT ( l2012h ; 
(17) IDessart et al.. (.2008) ; (18) Solanes et al. C2002) ; (19) Li et al.. (.2007) ; (20) .Maguire et al.. (.20im ; (21) ,Gil de Paz et al.. 0200% ! 

Notes. + From date of discovery (date of explosion is unknown). 
* Galactic reddening. 



dKrisciunas et al.', '2009') and fo = JD 2 452 859 as the moment 
of explosion (Jones et al., 2009), while we use E{B - V),otai 
- 0.19 mag as the estimated reddening of the host galaxy 
dOlivares et al l l2oToh . 

SN 2003ie was discovered in NGC 405 1 on 2003 September 
19.8 our & Boles, 200 3J) with an app arent brightness of ~ 
15.0 mag (un filtered). Benetti et all(l2003h classified it as a type- 
II SN, while iHarutvunvan et alF i 2008 ) found similarities be- 
tween some spectral features of SN 2003ie and the peculiar type 
II-P SNe 1998A and 1987A. Based on data on both the SN 
and i ts host galaxy (distance, metallicity, extinction) Smartt et al.l 
(12009.) gave 25 as an upper limit fo r its progenitor mass . 
For our stud i es we adopt D = 15.5 Mpc (iPierc e & Tull vlll988 : 
Smartt et al.l |2009b. E{B - V),^/ = .013 mag (Sc hlegel etaL , 
.1998), and fo = 2452 868 (IHarutvunvan et al.. 2008). 

SN 2004A was discovered in NGC 6207 on 2004 
Janu ary 9.8 by th e Japa nese amateur astronomer K. Itagaki 
(see iNakano et al.L |2004|) with a brightness of ~ 15.7 mag 
(unfilter ed). A preliminary optical spectrum obtained by 
[Kawakita et all (12004.) indicated that it was a type-II SN. Based 
on their detailed optica l photometric and spectroscopic analy- 
sis, iHendrv et aP (l2006l) found 2004A to be a classical type II- 
P SN with £) = 20.3 Mp c E{B - V),o,at = 0.06 mag, and to = 
2 453 01 1 . iMaguire et"al] (1201 Oh also published optical photo- 
metric and spectroscopic data from the photospheric phase. 

SN 2005ad was discovered in NGC 941 on 2005 
February 6.4 by K. I tagaki on unfiltered CCD-images (see 
iNakan o & Itagakit 12005) with a brig htness of ~ 17 . 4 mag . The 
very blue continuum on early spectra (iMorrell et al.L l20()5l) indi- 
cated that 2005ad was discovered in a very early phase, but the 
exact value of fo is still unknown. Lacking detailed studies, clas- 
sification of 2005ad is difficult, but some authors (Smar tt et al.L 
[2009; .Li et al... .2011.) re fer to it as a ty pe II-P SN. We use the 
value s ofD = 20.8 Mpc (iLi et al.Ll201 ll) and £(B - V)g„i = 0.06 
mag (ISchlegel et al.L ll998l) ! 

SN 2005af appeared in the nearby galaxy NGC 4945 on 
2005 February 8.2 with a n appa rent brightness of V = 12.8 
mag dJacques & PimentelL l2005l) . Spectroscopic analysis of 



iFiUppenko & Foley' ('2005') showed that it was a normal type II-P 
event. Pereyra et al. (2006) obtained optical polarir netry during 
the photos pheric and the early nebular phase, while iKotak et al.l 
(i2006i) and lKotakl(l2008l ) reported mid-IR data from this period 
(see details below). From the latter study we adopt D = 3.9 Mpc 
and fo = 2453 379, and use E{B - V)gai = 0.183 mag for redden- 
ing (Schlegel et al., 1998). 

SN 2005cs was discovered in the nearby, spectacular 
Whirlpool Galaxy (M51, NGC 5194) on 2005 June 28.9 by 
Kloehretal. (2005). It was classified as a type-II SN by 
iModiaz et al.l (l2005l) . Because of its proximity, it was the tar- 
get of many observing campaigns. Early-time optica l pho- 
tometry a nd spectroscopy wer e presented by Pastorell o et al.l 
(I2006h and lDessart et al.l (l2008h . while lBrown et al.l (2007) stud- 
ied the SN using the Swift ultraviolet-optical and X-ray tele- 
scopes. SN 2005cs seems to belong to the subtype of low- 
luminosity, ^^Ni-poor , low-velocity SNe II- P, as shown also by 
long -term studies of Tsvetkov et al.l (l2006h (optical photome- 
try), Gnedin et al. (2007) (R-band photometry and polarime- 
try), and .Pastorello et al.. (2009) (optical and near-IR photom- 
etry, optical spectroscopy). Extensive studies of the first-year 
data and pre-explosion images resulted in a progenitor mass of 
7-13 Mo (Maund et al., 2005; Li et al., 2006; Takats & Vinkl 
i2006l:IEldridge et aHI2007l) . while lUtrobin & Chugail ( 12008 ) de- 
termined a s ignificantly higher valu e (17-19 Mq). In their more 
recent work iPastorello et al.l (l2009l) concluded that the progen- 
itor mass should be in the range of 10-15 Mq. The distance to 
the host galaxy, deter mined via different methods, is scattered 
between 5.9-9.4 Mpc ("Feldmeieret al., 1997; Past orello et all 
20061: iTakats & Vinkd. i2006: .Dessart et al., .2008,) . Recently, 



Vinko et al I (120121) determined its mean value of 8.4 Mpc by 
combining the data of SNe 2005cs and 2011dh in the same 
galaxy. We adopt this latter value for D as well as fo = 
2453550fetorello et al., 2006), and E{B - V),o,„i = 0.04 mag 
(iBaron et aI.'U2007; .Dessart et al.,.2008,) . 

SN 2006bc was discovered in NGC 2397 by R. Martin (see 
iMonardetall l2006h on 2006 March 24.6; its apparent bright- 
ness was 16.0 mag (on an unfiltered CCD-image). Early spectra 
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clearly showed that it was a type-II SN (iPatat et all l2006t) . but 
due to the lack of long-term studies it remains unsolved whether 
it is a I I-P or a II-L SN (Smartt et al., 2009; Otsuka et al., 2Q\%. 
While I Smartt et al.l (l2009l) could not identify t he pr ogenitor of 



galaxy (iGil de Paz et all l2007h and E{B - V)gai 
(ISchlegel et allll998h as the galactic reddening. 



0.023 mag 



While lamartt et ai.l (izuu^) coum not laentiiy the progenitor or „„ , ., ^ . 

2006bc on HST archive images, lOtsuka et al.l (120121) als o used ^.2. Mid-infrared photometry with Spitzer 



HST to obtain late-time optical photometry. iBrown et al.l (l2009t) 
also observed it with the Swift UVOT detector. Parallel to our 
work, Gallagher et al. (2012) carried out a complete optical and 
infrared analysis (see details below). The distance to the host 
galaxy, D = 14.7 Mpc, is only estimated from its kinematic mo- 
tion (S martt et all l2009l). We also adopt E(B - V)^^/ = 0-205 
mag (ISchlegel et all 1 19981) . and fo = 2 453 8 19 dOallagher et all 
I2012h . 

SN 2006 bp, appearing in NGC 39 53, was found by K. 
Itagaki (see Nakano & Itagakil 12006a') on 2006 April, 9.6 
(16.7 mag , unfiltered). Early spectroscopic (Oui mbv et alll2006l 
1200 7*: Dessart et al., 2008) and p hotometr ic observations - in op- 



tical (.Ouimbv et al. 2007; Dessart et al.L I2OO8I). UV and X-r ay 



iby 

(llmmler et all 120071: iDessart et all l2008t iBrown et all l2009l) 
showed that 2006bp was a classical type II-P SN, and it was 
discovered only a few hours after shock breakout, which is 
the earli est discovery amon g type II-P SNe. Adopting the re- 
sults of Dessart et al.1 (l2008h . we used the data D - 17.5 Mpc, 
E{B - V),o,ai = 0.400 mag, and fo = 2 453 833. 

SN 2006my was discovered in NGC 465 1 by K. Itagaki (see 

iNakano & Itagaki, 2006b) on 2006 November 8.8 with an appar- 
ent unfiltered brightness of 15.3 mag, well past maximum light. 
Analyzing the first spectra after discovery, 'Stanishev & NielsenI 
(2006) classified it as a type II-P SN. Li et al. (200^ and 
[Maguire et alj (12010^ carried out optical photometric and spec- 
troscopic observations both in the photospheric and in the early 
nebular phase. From optical spectropolarimetrv. lChornock et alJ 
( I2OI O) showed strong late-time asphericities in the inner core 
of 2006my. Extended studies, based on the analysis of pre- and 
post-explosion images of HST and high-resolution ground-based 
telescopes, concluded that the progenitor of the SN sh ould have 
been a red supergiant with a mass between 7-13 (iLi et al 
20071: iLeonard et al', '20081 ISmartt et all 120091: ICrockett et al' 
201ll) .ISolanes et al. (2002) reported D = 22.3 Mpc as the dis 



tance of the host galaxy; we used th is value, as well as fn — 
2 453 943 as the moment of explosion dMaguire et all 201 Oh and 
E{B - V)gai - 0.027 mag as galactic reddening dSchlegel et all 
Il998h . 

SN 2OO60V, appearing in NGC 4303, has many simi- 
larities with 2006my. It was also found by K. Itagaki (see 
iNakano & Itagakil l2006cl) on 2006 November 24.9. Its unfil- 
tered apparent brightness was 14.9 mag at the moment of 
discovery. Similarly to 2006my, 2006ov was a type II-P SN 
disco v ered near the end of the plateau phase (iBlondin et al.L 
l2006t iLi etall l2007h. It was a lso one of the three CC SNe 
in which IChornock et al. (l201(?) found late-time asphericity in 
;t al.l T2007) determin ed the mass of the 
but later works dSmartt et all 120091: 
as an upper limit. We adopt 
and E{B-V)gai = 0.022 mag 



Li et i 



the inner ejecta. 
prog enitor as 15^, M0, 
ICrockett et al.„ . 201 1.) gave 10 Mq 

p = 12.6 Mpc (ISmartt et alll2009t) „ 

dSchlegel et all 1 1998*), while we estimated fo = 2 453 964 based 
s of |Li 



on the results of|Li et al. (2007). 

SN 2007OC was discovered in NGC 741 8 A on 2007 
Nove mber 3.1 by the Chile an Automatic Supernova Search pro- 
gram dPignata et all |2007|) . Due to the lack of any published 
photometric or spectroscopic observations, the moment of ex- 
plosion and its subtype remain unknown; it is me ntioned in the 
ASIA GO Supernova Catalogue as a type II-P SN dBarbon et all 
l2008l) . We adopted D - 28.0 Mpc as the distance of the host 



We collected and analyzed all available IRAC post-basic cali- 
brated data (PBCD) on the 12 CC SNe Usted above. The scale 
of those images is 0.6"/pixel. In some cases, a clear identifica- 
tion of SN was not easy because of the faintness of the target 
in the vicinity of bright regions. We tried to apply the image- 
subtraction technique us ed succesfully in some cases t o ana- 
lyze S pitzer data of SNe dMeikle et all[200 6. 2007; Kota k et all 
l2009i) . We found pre-explosion images of the host galaxies for 
only two SNe (2005cs and 2006ov), and their quality was too 
poor for reliable image subtraction. We also generated some sub- 
tracted images using the latest frames as templates, but their ap- 
plicability was also very limited. Finally, we were able to un- 
ambiguously identify a MIR point source in nine cases at the 
SN position, while in the cases of three SNe (2003hn, 2005cs, 
and 2006bc) we did not find any sign of a point source, either 
on the original images or on the subtracted image (however, as 
we mentioned before, there are some other recent results about 
SN 2006bc, see details in Section|4|. For the nine detected SNe 
(Fig-[D we carried out a detailed photometric analysis. 

Originally, we planned to apply PSF-photometry using the 
DAOPHOT package of IRAlQ to obtain more accurate results. 
Unfortunately, the bright background flux of the host galaxies, 
sometimes showing rapid spatial variations, and the lack of rel- 
atively bright point sources on the frames prevented us from ap- 
plying this. 

Therefore we carried out simple aperture photometry on the 
PBCD frames with the PHOT task of the IRAF software pack- 
age. This was possible beacuse the studied SNe - except of 
the three mentioned above - were clearly identifiable on IRAC 
images (see Fig. [1]). We generally used an aperture radius of 
3.6"and a background annulus from 3.6"to 8.4"(3-3-7 configu- 
ration in native ~ 1 .2"IRAC pixels), applied aperture corrections 
of 1.124, 1.127, 1.143, and 1.234 for the four IRAC channels 
(3.6 , 4.5, 5.8 , and 8 .0 jum) as given in the IRAC Data Handbook 
(Re ach et all l2006h . We chose a 2-2-6 configuration for SNe 
2003J and 2007OC (aperture corrections: 1.213, 1.234, 1.379, and 
1.584), and a 5-12-20 configuration on the images of the bright, 
nearby SN2005af (aperture corrections: 1.049, 1.050, 1.058, and 
1.068). 

There is one object, 2005af, f or which MIR data were partly 
analyzed and published before: iKotak etakl (|2006|) presented 
IRS spe ctra and also IRAC and PUI fluxes up to -1-433 days, 
while in [Kot^(f2008b there is a figure showing an IRS spectrum 
and IRAC measurements carried out at -1-571 days. To check 
the reliability of our reduction technique we also re duced these 
data and found a good agreement with the results of Kot ak et al] 
(2006) and Kotak (2008) . We also completed this dataset with 
the fluxes determined from the latest measurements found in the 
Spitzer database. 

We also carried out aperture photometry using IRAF on late- 
time MIPS 24.0 yum PBCD frames (image scale 2.45 "/pixel) for 
SNe 2003J, 2003ie, 2005af, 2006bp, and 2006my. We used an 
aperture radius of 3.5"and a background annulus from 6"to 8", 
except for SN 2005af, for which we found the 5-5-12 configura- 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 2. IRS spectra of SN 2005 af. The spectra at different epochs 
are vertically shifted relative to each other by the amount indi- 
cated in the labels next to the spectra. 



tion to be the best choice. For SNe 2005af and 2006my, the MIR 
point sources were also identifiable on available IRS PUI im- 
ages. Their fluxes were measured using the MOsaicker and Point 
source Extractor software. All measured fluxes are 

listed in Table lAl] 

2.3. Mid-infrared spectroscopy witli IRS 

We also downloaded the available IRS spectra of the SNe from 
the Spitzer database. Four spectr a of SN 2005af - th e firs t three 
have a lready been published by iKotak et aP (l2006t) and iKotakI 
(120081) -, two spectra of SN 2004A, and a single spectrum of 
SNe 2003ie, 2005ad, 2005cs, 2006bp, and 2006ov were col- 
lected this way. 

We analyzed the PBCD-frames of the spectroscopic data 
with the SPitzer IRS Custom Extraction software {SPIC^. Sky 
subtraction and bad pixel removal were performed using two ex- 
posures that contained the spectrum at diff'erent locations, which 
were then subtracted them from each other. Order extraction, 
wavelength- and flux-calibration were performed by applying 
built-in templates within SPICE. Finally, the spectra from the 
first, second and third orders were combined into a single spec- 
trum with the overlapping edges averaged. In a few cases the sky 
near the order edges was oversubtracted because of some excess 
fluxes close to the source, which resulted in spurious negative 
flux values in the extracted spectrum. These were filtered out us- 
ing the fluxes in the same wavelength region that were extracted 
from the adjacent orders. 

We found that the spectra of SNe 2005cs and 2006ov are 
too noisy for further analysis. For the other SNe the spectral and 
photometric fluxes agree reasonably well. We also compared the 
extracted spectra with the public data in the Cornell Atlas of 



^ http://ssc.spitzer.caltech.edu/dataanalysistools/tools/mopex/ 
^ http://ssc.spitzer.caltech.edu/dataanalysistools/tools/spice/ 



Spitzer IRS Sources Q (CASSIS, iLebouteiller et al.L l201lh . and 
found them to be similar to each other. 

The spectra of SN 2005af are shown in Figure El The three 
ear hest spectra w ere already published by iKotak etal] (l2006h 
and iKotakI (|2008|) . Our spectra agree well with theirs, which 
serves as a good test for checking our data reduction. In Figure 
[3] we present the extracted IRS spectra for the other SNe. 
Although detailed modeling of these spectra is beyond the scope 
of this paper, we identified spectral features of SiO, and the 

I- vibrationa l transi tion of CO based on fi ne identifications 
of IKotak et all (l2005h and lKotak etaH (l2006l) . The SiO feature 
plays an important role in modeling the observed SEDs, while 
the CO 1-0 transition results in a flux excess at 4.5 micron (see 
Section[3]l. 

3. Analysis 

3.1. Mid-IRSEDs 

Figure |4] shows the mid-IR SEDs of the SNe calculated from 
IRAC 3.6, 4.5, 5.8, and 8.0 jum fluxes (plotted with available PUI 
13.0-18.5 fim and MIPS 24.0 //m fluxes). The continuum fluxes 
of the reduced IRS spectra match the SED flux levels well. 

As is seen in Fig. H IRAC fluxes of several SNe (2003J, 
2004A, 2005ad, 2005af, 2007oc) show obvious variations in 
time, while others do not. Based on t he a nalysis of other 
CC SNe, e. g., SN 2004et (Kotaket SI 120091) or SN 2004dj 
(ISzalai et allboi UlMeikle et aU,201 i fTthis temporal evolution 
could be the sign of newly-formed dust grains in the ejecta, but 
there may be more effects that have to be taken into account to 
explain this phenomenon. 

3.2. Fitting of dust models to MIR SEDs 

Assuming that the radiation is purely thermal, the main source 
of mid-IR flux is most likely warm dust. To estimate the min- 
imum amount of dust, we fit analytic dust models to the ob- 
served SEDs. Prior to fitting, the observed SED fluxes were 
dere ddened using th e galac tic reddening law parametrized by 
Fitz patrick & Massal (l2007l) assuming Ry - 3.1 and adopting 
E(B - V) values as listed in Section lTTI 

As a first step we examined whether simple, pure blackbody 
(BE) emission might be consistent with the observed MIR SEDs. 
Although a pure BE model has little physical significance, we 
applied it to estimate the minimum radius of an optically thick 
dust sphere that fits the observed SEDs. From the minimum 
radii we calculated the corresponding velocities assuming ho- 
mologous expansion (see Table[3]l. The typical velocities of type 

II- P SNe are 2-3000 km s ' for the inner ejecta where dust for- 
mation may take place. Therefore, for a larger minimum size it 
is not possible that the inner ejecta reach the calculated size of 
the EE. This means that in these cases the minimum radii of the 
optically thick dust spheres are too large and that consequently, 
most of the estimated dust mass cannot have originated from 
newly formed grains. 

The analytic models were calculated using Eq.l in 
Meikle et al. (20.(11) assuming a homogeneous (constant- 
density) dust distribution: 

= 2t?R^B,(T) [t;2(2t2 - 1 + (2t, + 1) exp(-2T,)] , (1) 

The Cornell Atlas of Spitzer/IRS Sources (CASSIS) is a product of 
the Infrared Science Center at Cornell University, supported by NASA 
and JPL. 
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Fig. 1. Type II-P SNe on IRAC post-BCD 4.5 jum images. The FOV is 100"x 60". North is up, and east is to the left. 
Table 2. IRS observations of the SNe 



Object 


UT Date 


MJD - 




ID 


Proposal ID (PID) 






2450000 


(days) 






SN 2003ie 


2005-01-14 


3385 


517 


rl0557696 


3248 (Meikle et al.) 


SN 2004A 


2004-08-27 


3244 


233 


rl0557952 


3248 (Meikle et al.) 




2005-02-16 


3417 


406 


rl0558208 


3248 (Meikle et al.) 


SN 2005ad 


2005-09-08 


3621 


215+ 


rl4466304 


20256 (Meikle et al.) 


SN 2005af 


2005-03-17 


3446 


67 


rl3413376 


237 (Van Dyk et al.) 




2005-08-11 


3593 


214 


rl4468096 


20256 (Meikle et al.) 




2006-08-03 


3950 


571 


r 17969664 


30292 (Meikle et al.) 




2007-03-25 


4185 


806 


r 17969920 


30292 (Meikle et al.) 


SN 2006bp 


2008-01-18 


4483 


648 


r231 11936 


40619 (Kotaketal.) 



Notes. From date of discovery. 



where R is the radius of the dust sphere, By{T) is the Planck 
function, T is the temperature of the dust, and Ty is the op- 
tical depth at a given frequency. To estimate the dust opti- 
cal depth, we adopted the p o wer-la w grain-size distribution of 
iMathis. RumpI & Nordsieckl (Il977h . hereafter MRN, assuming 
m = 3.5 for the power- law index and grain sizes between a„„„ - 
0.005 fim and a„,ax = 0.05 ;um. Ty can then be calculated as 



1 



(2) 



where pgrain is the density of dust grain material, k is the grain 
number-density scaling factor, and Ky is the mass absorption co- 
efficient. The dust is assumed to be distributed uniformly within 
a sphere. R, T, and k were free parameters during the fitting. 



First we calculated the dust models using amorphous 
carbon (AC) grains. Dust opacity values were taken from 
Colangeli et al. ( 1995), while the density of t he dust grain mate- 
rial, Pg rain - 1-85 g cm"^, was adopted from lRouleau & Martini 
(Il99lh . Although we were able to successfully apply these 
models in most cases, finding adequate solutions for SNe 
2005af and 2006my was possible only by assuming silicate dust 
(C-S i-PAH mixture, MRN-d istribution, pgrain = 1-85 g cm"^, 
We ingartner & DraineL I2OOII) . The IRS sp ectra of SN 2005a f 
also show the signs of SiO (see Sectionl2land lKotak et al.Ll2006l) : 
but there were no IRS spectra on SN 2006my in the database. 
On the other hand, the spectra of the other SNe (2003ie, 2004A, 
2005ad, 2006bp) that can be fit well with AC models do not 
show obvious signs of SiO. The parameters of the best-fit ana- 
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Fig. 3. IRS spectra of SNe 2003ie, 2004A, 2005ad, and 2006bp. The spectra of SN 2004A are vertically shifted relative to each other 
by the amount indicated in the label next to the upper spectrum. Strong and wide emission features in the spectrum of SN 2006bp 
possibly emerge from nearby interstellar clouds. 



lytic models (shown in Figure |5]l are collected in Tables |4] and 

12 

In a few cases we found that a single temperature model can- 
no t describe the whole o bserved SED. This is similar to the result 
of IWooden et~an (Il993h . who showed that for SN 1987A a hot 
component (r ~ 5- 10 000 K) may affect the continuum emission 
of the warm dust. This hot component is thought to be caused by 
an optically thick gas in the innermost part of the ejecta. The con- 
tribution of such a hot component can be significant at shorter 
IRAC wavelengths. The modeling of the hot component needs 
late-time optical or near-IR data taken simultaneously with the 
Spitzer observations. Because of the lack of such data for most 
of the SNe in our sample, we usually did not include the 3.6 ^m. 
(and sometimes also 4.5 /im) fluxes while fitting the dust models. 
There are two exceptions for which we have found optical pho- 
tometry obtained close to the IRAC measurements: SN 2004A 
at 247 days (Hendry et al.L [2006), and SN 2006my at 205 days 
jMaguire et akl 20101) . At these epochs we completed the ob- 
served SEDs with the optical data and added a hot component 
to our models (see the fittings in Figure|6] and the parameters in 
Table|6]). 

Figure |5] shows that SNe 2004A, 2005ad, 2005af, and 
2007oc show variable excess flux at 4.5 jim (in these cases we 
omitted the 4.5 jim fluxes from the fitting). The excess flux 
weakens with age, as the ejecta temperature decreases. This 
tren d is also seen in the IRS spectra of SN 2005af (see Section|2] 
and iKotaket al.L 12006). This may be due to the influence of the 
red wing of the 1-0 vibrationa l transition of CO at 4.65 jjm (see 
also iKotak et all I2005L 120061: ISzalai et all 1201 Ih . An alterna- 
tive explanation for temporarily changing flux at 4.5 may 



be the rapid destruction of stochastically heated grains. This ef- 
fect becomes_si^ificant when the grain sizes are below 0.003 
lim dBocchio et al.L l20T2h . However, as we mentioned above, in 
the cases studied here the grain sizes are thought to be more 
than 0.005 j um (s ee also Section [1] and Kozasa et all 120091 : 
I Szalai etall 1201 Ih . Therefore the rapid grain destruction pro- 
cess seems to be less significant around these young SNe than it 
might be in older SNRs. 

We also found that single-temperature models underestimate 
the fluxes at 24 /vm when there are MIPS data taken contempo- 
raneously with the IRAC measurements (SNe 2003ie, 2006bp, 
2006my). In these cases w e added a co l d com pon ent {T ~100 
K) to o ur models, similar to Kotak et alj ( l2009l) and lSzalai et al.l 
which resulted in a reasonably good fit to the whole ob- 
served SED. The parameters of the cold EEs are collected in 
Tabled 

The calculated dust masses for all SNe and all epochs are 
also collected in Table |4] and |5] Masses were determined using 
the formula 

Md = AnR\yl2,K^ (3) 

adopted from iLucv et alJ ([1989') and iMeikle et al.1 (l2007h . It 
should be no ted that thes e masses from analytic models are 
lower limits (iMeikle et all 12007 ). because of the implicit as- 
sumption that the dust cloud has the lowest possible optical 
depth. As mentioned before, more realistic dust masses can be 
obtained only from numerical models for the optically thick dust 
cloud. 

To check this statement and the trustworthiness of our 
analytic model fitting, we generated some models with the 
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Fig. 4. Mid-IR SEDs of the nine CC SNe at different epochs. Open symbols denote IRAC warm mission data. 



three-dimensional radiative-transfer code MOCASSIN 2.02.55. 
(Ercolano et al., 2003, 2005). This numerical code uses a ray- 
tracing technique, following the paths of photons emitted from 
a given source through a spherical shell containing a specified 
medium. The shell is mapped onto a Cartesian grid allowing 
light-matter interactions (absorption, re-emission, and scattering 
events) and to track energy packets until they leave the shell and 
contribute to the observed SED. It is applicable for reconstruct- 
ing dust-enriched environments of CC SNe and to determine 
physical parameters suc h as grain-size di s tribution, composition , 
and mass of dust (see Sugerman et al., 2006!; lErcolano et al.L 
[2007; Fabbri et al.L l20il; Szalai et al., 201 1). The code handles 
different types of dust heating processes caused not only by pho- 
tons but also by collisions of grains and electrons (B. Ercolano, 
personal comm.). 

We chose SN 2005af as a test object (at 576 days), for which 
we could fit both AC and Si-containing dust models. For AC 
grain composi tion we used the built-in optical constants from 
iHannej (11988 ). while to fit also the flux at 8 fim we found a 
0.65:0.35 Si-AC mixture to be the best solution (built-in optical 
cons tants of astronomical silicate are taken from lLaor & DraineL 
Il993i) . We applied the MRN grain-size distribution. The final 
results of numerical calculations are shown in Table |7l while 
we compare the best-fitting analytical and numerical models in 
Figure|2] The excess flux below 3 jum in the MOCASSIN model 
comes from the assumed hot central source, which was absent in 
the analytic models. 



Although the two methods are very different (the analytic 
models describe the thermal MIR radiation of the dust, while 
the calculations of MOCASSIN are based on more detailed ra- 
diative transfer), we were able to generate similar SEDs in the 
MIR region with similar dust shell sizes. For pure AC dust, dust 
mass derived with the numerical model is an order of magnitude 
higher than that derived with the analytical model, which is ver y 
similar to our earlier results on SN 2004dj (ISzalai et al.Ll201 11) . 
Models for silicate dust are slightly different in numerical and 
in analytic models, and, to obtain the best fit, we had to apply 
larger grains (a = 0.1-1 fim, MRN distribution). Despite these 
differences, the calculated mass for the silicate dust shell is very 
close to the value given by the analytic model. 

4. Discussion 

To understand the origin and heating mechanisms of dust in the 
environment of SNe, it is important to combine t he MIR data 
with optical and/or near-IR measurements (see e. g. Kotak et aU 
2009; Foxetal., 2010; OtsukaetaL, 2010; Szalai et al., 201ll 
iMeikle et al.L 1201 L) . Observations in X-ray and radio may also 
provide useful additional information about SN-CSM interac- 
tions. However, while Spitzer observations on the studied SNe 
were typically carried out in the nebular phase, there are only 
a very few additional late-time observations of these objects 
(see these cases above). Therefore, we can draw our conclusions 
based on mainly the MIR data analyzed in the previous Section. 
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Fig. 5. Best-fit analytic models on MIR SEDs of the SNe. The 3.6 jum and sometimes the 4.5 //m fluxes are usually not included in 
the fitting (see text). The fittings that contain hot and warm components for SN 2004A and SN 2006my are shown in Figure|6] 
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Fig. 6. Best-fit models for SN 2004A at 247 days (left) and for SN 2006my at 205 days (right) that contain a warm dust component 
and a hot blackbody (see text). 



As mentioned before, there are three SNe in our sample 
(2003hn, 2005CS, and 2006bc) for which the available MIR data 
were not appropriate to resolve the point source at the position 
of the SN. The case of SN 2006bc is special, because there is 
a bright source (probably a surrounding H II region) very close 
to the position of the SN. One would need high-resolution data 
and more sophisticated methods to extract inf ormation on this 
SN (lOtsukaetal .1120121) . Parallel to our studies. [Gallagher et all 



(120 12h carried out a complete optical and infrared data analysis 
of SN 2006bc based not only on Spitzer observations but also 
HST and Gemini South data. They found M^,,,, < 2 x 10"^ Mq 
newly formed dust in the ejecta, which is consistent with other 
similar cases. 

The SEDs of the studied SNe and the derived parameters 
shown in Tables |3]|4] and|5]give guidelines and limits to the pos- 
sible sources of the detected MIR radiation. We have the highest 
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Table 7. Parameters for the best-fit MOCASSIN models to SN 2005af SED at +576 days. 



Composition 






^our 






of dust 


(10^ Lo) 


(K) (10'«cm) 


(10"* cm) 


10-« cm-3 


(10-5 Mo) 


AC 


3.4 


7000 


0.8 


150+ 


18.0 


SiAC (0.65:0.35) 


16.9 


8000 0.2 


1.0 


0.02* 


11.0 



Notes. T, and L, denote the temperature and the luminosity of the central illuminating blackbody source; the geometry of the dust cloud was 
assumed to be a spherical shell with inner and outer radii i?,„ and i?„„,. 
t Homogeneous (constant-density) spatial distribution of grains 
* p oc r-' density profile for spatial distribution of grains 



Table 3. Minimum radii and corresponding velocities deter- 
mined from the best-fit BB-models to observed SEDs. 



Object 


Epoch 




Vbb 




(days) 


(10'' cm) 


km s ' 


SN 2003J 


47r 


6.3 


15 480 


SN 2003ie 


488 


6.0 


14 230 




632 


4.9 


8975 




1763 


6.1 


4005 


SN 2004A 


247 


3.4 


15 930 




445 


3.8 


9885 




563 


2.2 


4520 


SN 2005ad 


198' 


0.3 


1755 




364' 


0.2 


635 


SN 2005af 


194 


0.8 


4770 




399 


0.4 


1160 




576 


0.6 


1205 




772 


0.8 


1200 




940 


0.6 


740 


SN 2006bp 


401 


4.5 


12 990 




628 


5.5 


10140 




767 


5.1 


7695 


SN 2006my 


205 


3.8 


21455 




342 


3.7 


12 520 




559 


4.2 


8695 


SN 2OO60V 


321 


5.1 


18 390 




539 


3.4 


7300 


SN 2007OC 


250' 


5.7 


26285 




415' 


4.9 


13 665 



Notes. * From date of discovery. 



amount of data for SN 2005 af, for which a sequence of MIR 
SEDs is available that shows a monotonically decreasing tem- 
perature (see Figure |4]i. This may be the sign of newly formed, 
cooling dust. It is also supported by the size of the dust sphere 
determined from the models: the calculated expansion velocities 
are lower than 1300 km s-' in the nebular phase, which is low 
enough to be compatible with the inner ejecta velocities of II- 
P SNe. This is clearly another argument that most of the dust 
may have been formed in the ejecta. The calculated dust masses 
(1-3 X 10-"* M0) agree with the amount of SiO calculated by 









SN 2005af +576 days 








450KC-Si-PAH 

Mocassin 


CVJ 


1e-17 






'e 










(erg s 








Flux 


1e-18 







3 5 10 20 

Wavelength (\x.m) 



Fig. 7. Best-fitting Si-containing analytical and numerical mod- 
els for SN 2005af at -1-576 days (solid and dotted lines, respec- 
tively). The excess flux below 3 ;um in the MOCASSIN model 
is due to the assumed hot central source, which was absent from 
the analytic models. 

iKotak et al.1 (l2006h and wifli the results of lKotakI (l2008h . and they 
are also similar to the values found in other SNe. 

Another explanation for the observed behavior of MIR SEDs 
may be the thermal radiation of pre-exis ting dust in the CSM 
re-heated by the S N, i. e., t he IR-echo ( Bode&Evansl 1 19801: 
Ibwek, 1983, 198^ ISugerma n. 2003). We have very httle in- 
formation on the properties of the CSM around this SN. There 
are no published X-ray or radio measurements, and the rela- 
tively high value of the interstellar reddening (E(B — V)gai = 
0.183 mag) is caused by the ga lactic ISM. However, based on 
iDwekl ([1983) and Meikle et all (i2006), the shapes of the MIR 
light curves should be rather flat if there is an IR echo, but 
the IRAC light curves of SN 2005af decline monotonically in 
time (see Figure [8]l. Therefore we conclude that an IR echo is a 
less probable source of MIR radiation during the observed inter- 
val. It is also possible that dust grains condense in a cool dense 
shell (CDS) that is generated between the forward and reverse 
shock waves during the interaction of the SN ejecta and the pre- 
existing CSM. This is usually invoked for type Iln SNe, but it 
was also taken into account in a few other cases to explain the 
presence of early or very late MIR-exce ss dKotak etaU 120091: 
[Andrews et al 1 12OIOI: iMeikle et aI.L l201 ll) . A CDS was also as- 
sumed to be the source of the cold dust component in the SEDs 
of II-P SNe (Szalai et al., 2011). 

The shapes of the MIR light curves and the presence of SiO 
make SN 2005af very similar to the well-studied SN 2004et. In 
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Table 4. Parameters for the best-fit analytic amorphous carbon 
dust models to the warm components of SEDs. 



Epoch 7^vvc7/7» ^wanu ^ ^dust ^warm 

(days) (K) (10'" cm) (IQ-^^) (IQ-^ Mq) (lO^'ergs-') 

SN 2003J 

471'' 370 6.4 54.5 710 48.8 

SN 2003ie 

488* 310 6.2 60 700 22.2 
632 300 10.1 4 200 16.9 
1763 280 15.8 L5 290 19.1 

SN 2004A 

247* 340 3.4 100 160 8.5 
445 310 3.9 67 200 8.1 
563 370 2.2 350 180 6.2 

SN 2005ad 

198* 890 0.4 480 1.0 4.9 
364* 750 0.4 190 04 1.8 

SN 2005af 

194 590 2.3 7.0 4.2 10.2 

399 640 0.6 190 1.6 3.3 

576 460 0.9 85 2.6 1.5 

772 430 0.8 66 0.7 0.3 

940 380 0.7 260 4^5 0.6 

SN 2006bp 

401* 370 4.6 100 480 26.2 
628 330 6.1 89 1000 29.3 
767* 350 5.2 100 690 27.0 

SN 2006my 

205* 470 3.8 110 88 19.3 
342* 420 4.1 67 230 32.5 
559* 410 4.3 93 370 34.3 

SN 2OO60V 

321 350 5.1 100 650 25.9 
539 350 8.2 31) 82 15.4 

SN 2007OC 

250* 340 8.9 9.0 310 37.8 

415* 340 5.5 45 370 23.5 



Notes. * From date of discovery. 

* An additional (cold or hot) component was also fit (see details in Table 
[6] and in the text). 



this SN the main source of the warm component is also iden- 
tified as coohng, newly formed dust in the ejecta (iKotak et al.L 
|2009 |). Unfortunately, there are no Spitzer measurements on 
SN 2005af at very late epochs, so we could not see whether there 
is also a rise in IRAC fluxes after -hIOO O days like in SN 2004et 
(IKotak et all 120091: iFabbri etaLl l20Tll) . 

There are fewer Spitzer data for the other eight SNe, and 
some of them were recorded during the warm mission (since the 
end of the cryogenic mission in 2009 only the IRAC 3.6 and 4.5 
/im channels have been usable). Nevertheless, we may also draw 
some conclusions for these SNe. Among them we found obvious 



Table 5. Parameters for the best-fit analytic silicate dust models 
to the warm components of SEDs. 



Epoch ^wflrm ^warm ^ ^dusj ^warm 

(days) (K) (10'* cm) (lO^-^) (10"^ Mq) (lO^'ergs-') 



SN 2005af 

399 550 1.1 200 21 4.3 

576 450 1.0 450 31 2.0 

772 390 0.8 410 17 0.7 

940 400 0.8 160 8J a5_ 

SN 2006my 

205" 380 8.0 22 920 44.9 

342"'" 400 5.5 100 1300 33.8 

559"^ 380 6.2 100 1900 28.4 



Notes. ' An additional (cold or hot) component was also fit (see details 
in Table|6]and in the text). 



Table 6. Parameters of the additional cold and hot blackbodies. 



Object 


Epoch 


T 


R 


L 




(days) 


(K) 


(10'* cm) 


(10^' erg s-') 


Cold component 


SN 2003ie 


488 


90 


29.5 


11.6 


SN 2006bp 


401 


110 


24.0 


6.0 




767 


80 


42.0 


5.1 


SN 2006my 


342 


120 


33.5 


18.0 




559 


90 


76.0 


28.6 


Hot component 


SN 2004A 


247 


4800 


0.045 


76.1 


SN 2006my 


205 


3700 


0.105 


146.3 



MIR flux changes in four objects: SNe 2003J, 2004A, 2005ad, 
and 2007oc (we note that for SN 2003J there is only one com- 
plete IRAC measurement, the second one was made during the 
warm mission). From these, SN 2005ad is the only object for 
which calculated parameters may indicate dust formation with 
high probability, because between the two observed epochs the 
dust model temperature decreased from 890 K to 750 K. The cal- 
culated expansion velocities are even lower than for SN 2005af, 
and the calculated dust masses (~10'^ M0) also agree with the 
values found in similar cases. Interestingly, for both SNe 2005af 
and 2005 ad, the calculated radii of the dust sphere are roughly 
constant (or decline shghtly) in t he nebular phase in spite of the 
expected expansion. iKotak et al.l (12009) found a similar effect in 
SN 2004et: their explanation was that the dust was contained 
within an optically thin cloud of optically thick dust clumps that 
cannot move outward because of the hotter, higher pressure in- 
terclump gas. 

While the flux changes of SNe 2003J, 2004A, and 2007oc 
suggest that the observed MIR radiation is connected to the SN 
(at least partly), only a part of it may have been emitted by newly 
formed grains in the ejecta. The derived dust sphere radii indi- 
cate very high expansion velocities that are much higher than the 
inner ejecta velocities of type II-P SNe in the nebular phase. At 
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Table 8. List of type II-P supernovae with published Spitzer data. SNe typed with italic were analyzed by our group (this is only 
pai-tly true for SNe 2004dj and 2005af). 



Name 


Spitzer data 


Epochs 


Detectable local dust 


Dust mass References 










(10-= Mo) 



SN 2002hh 


IRAC, MIPS 


590,684,758 


Yes - IR echo 


< 3600 


1, 2 


SN 2003gd 


IRAC 


496-499,670-681,1264 


Yes - new, IR echo 


10 - 2000 


3,4 


SN 2004dj 


IRAC, MIPS, IRS 


98- 1381 


Yes - new, CDS, (maybe IR echo) 


1-80 


5, 6 


SN 2004et 


IRAC, MIPS, IRS 


65 - 1240 


Yes - new, IR echo 


4-200 


7, 8 


SNlOOSacT 


IRAC, IRS 


198, 215, 364 


Yes - new 


0.5-1 


tliis work 


SN 2005af 


IRAC, MIPS, IRS 


194 - 940 


Yes - new 


8-21 


9, 10, this work 


SN 2006bc 


IRAC, MIPS 


537-540, 678-687, 753-775 


Yes - new 


< 200 


11 


SN 2007it 


IRAC 


351,561,718, 944 


Yes - new 


5-70 


12 


SN 2007od 


IRAC, MIPS, IRS 


300, 455-470, 667 


Yes - new (maybe CDS) 


17-19 


13 


5^ 2003 f 


IRAC, MIPS 


471, 852, 2562 


Yes, partly 


< 710 


tliis work 


SN2004A 


IRAC, IRS 


233-247, 406, 445, 563 


Yes, partly 


<200 


tliis work 


SN 2007oc' 


IRAC 


251,415,634, 759 


Yes, partly 


< 370 


this work 


SN2003ie 


IRAC, MIPS, IRS 


474-488, 517, 632, 1763 


No 




tliis work 


SN2006bp 


IRAC, MIPS, IRS 


401-421, 628-671, 


No 




tliis work 






768-771, 1361 








SN2006my 


IRAC, MIPS 


205, 342-354, 559-568 


No 




this work 


SN2006OV 


IRAC, IRS 


321-343, 539 


No 




this work 



Re ferences. (l)lBarlow et al.l(l2005h: (2 )lMeikle etalj (120061): (3)[Su german et aU (120061): (4)jMeikle et al.U 2007 l) : (5) ISzalai et all (1201 ih: 
(6)|Meikle etal.'(2011); (7) Kotak et al. (2009): (8)' Fabbri et alj ( 1201 II) : (9) lKotak et alj J2OO6I) : (lO) lKotakj 120081) : (1 D iGallagher et alj ( 120121) : 



(12) Andrews et al. (2011); (13) An drews etaL(.2O10ir 
Notes. ' From date of discovery. 
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Fig. 8. IRAC hght curves of SN 2005af: 3.6 yum (open squares), 
4.5 /im (filled squares), 5.4 yum (open circles), 8.0 yum (filled cir- 
cles). 



such large radii grain condensation is unlikely. SN 2004A shows 
changes in the dust temperature, and the calculated luminosities 
of the dust spheres are not much higher than for SNe 2005af and 
2005ad, but the source of most of the MIR radiation may be pre- 
existing dust or the bulk of grains formed in a CDS. This may be 
even more valid for SNe 2003J and 2007 oc. 



For the remaining four SNe (2003ie, 2006bp, 2006my, 
2OO60V) there is no sign of real changes in the detected MIR 
fluxes. However, model fits to their SEDs suggest very high ex- 
pansion velocities and much larger dust masses than for SNe 
2005af and 2005ad (note that these values are still lower than the 
theoretically predicted dust masses). While this does not com- 
pletely rule out that there is newly formed or pre-existing dust in 
the close environment of these SNe, the amount of dust should 
be smaller than the calculated values. We conclude that there are 
probably other sources of MIR radiation, e. g., dust and/or cold 
gas clouds within the host galaxies or in the line-of-sight of these 
SNe. For SN 2006bp, the high value of reddening, E{B-V),o,ai = 
0.4 mag, supports this idea. For the other three SNe the E(B - V) 
values are much l ower but are es t imated only from the galactic 
reddening map of ISchle^el et all (!l998l). We also note that the 
3" aperture radii used by IRAC photometry correspond to 20-35 
pc at the distances of our SNe, which is 2-3 orders of magnitude 
larger than the regions that the SN ejecta or shock waves could 
achieve at the given epochs (it may be less valid for SN 2005af, 
which is much closer than the other SNe in the sample). 

In Table [8] we show all type II-P SNe that have Spitzer data 
analyzed and published to date including the objects presented 
in this work. As the final list shows, we doubled the number of 
objects, which might serve as a basis for statistical studies on a 
larger sample. We separated SNe that have some amount of local 
dust from those where there is no unambiguous sign for dust. 
Similar tables about the dust formation statistics of SNe can be 
found for instance in lFox et ah (i2011i) and iGall et ah (,201 L) . 

5. Summary 

We collected the public photometric and spectroscopic data of 
the Spitzer Space Telescope on 12 type II-P SNe. We found the 
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data of nine of these objects to be appropriate for a detailed 
study, which almost doubles the number of type II-P SNe with 
a detailed, published MIR data analysis. We fit the observed 
SEDs with simplified one- or two-component dust models. In 
SNe 2005 ad and 2005 af we found cooling temperatures and de- 
creasing luminosities of the warm component that are similar to 
the values found in other SNe that are thought to have newly 
formed dust in their environment. 

The calculated temperatures for the other SNe do not show 
strong temporal variation, while the derived luminosities and 
radii are too high to be compatible with local dust. Moreover, 
the calculated dust masses in these cases are orders of magnitude 
higher than the observed amount of dust around the well-studied 
SNe listed above. The large radius of the warm component may 
suggest pre-existing dust in these cases, making it unclear if 
there was new dust formed around these SNe. Nevertheless, the- 
oretical models predict orders of magnitude more newly formed 
dust in CC SNe. Our conclusions support the previous obser- 
vational results that warm new dust in the environment of SNe 
contributes only sUghtly to cosmic dust content. 

At the same time, our study may help to acquire better statis- 
tics on the dust formation capabilities of type II-P SNe, and 
present some objects that may be interesting for additional obser- 
vations at longer wavelengths. In a few years the SNe for which 
we concluded the presence of some local dust (2003J, 2004A, 
2005ad, 2005af, 2007oc) may show some other signs of dust 
formation in the transition phase. 
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T. Szalai and J. Vinko: Twelve type II-P supemovae seen with the eyes of Spitzer 



Table A.l. Spitzer photometry for the clearly identifiable supemovae in our sample. 











Flux (10-2" erg 


s ' cm 2 A 


') 


UT Date 


MJD - 






IRAC 




IRS PUI MIPS 




2450000 


(days) 


3.6 jum 


4.5 


5.8 


8.0 jum 


13.0 - 18.5 urn 24 urn 


SN 2003J 


2004-04-27" 


3123 


471+ 


1069(126) 


706(67) 


1311(95) 


1684(121) 




2005-05-13" 


3503 


852+ 










366(59) 


2010-01-17* 


5214 


2562+ 


750(116)) 


398(56) 








SN 2003ie 


2004-12-03^ 


3342 


474 










168(48) 


2004-12-17"' 


3356 


488 


152(44) 


77(23) 


363(36) 


515(38) 




2005-05-10"' 


3500 


632 


188(44) 


97(23) 


346(34) 


514(38) 




2008-06-14^ 


4631 


1763 


172(44) 


74(22) 


317(36) 


543(45) 




SN 2004A 


2004-09-10^ 


3258 


247 


287(48) 


257(34) 


127(23) 


162(19) 




2005-03-27"' 


3456 


445 


122(34) 


68(19) 


70(19) 


123(18) 




2005-07-23'' 


3574 


563 


93(30) 


47(16) 


78(19) 


110(19) 




SN 2005ad 


2005-08-23^ 


3605 


198+ 


200(40) 


149(26) 


106(17) 


47(9) 




2006-02-05" 


3771 


364+ 


71(26) 


81(20) 


43(12) 


20(7) 




SN 2005af 


2005-07-22" 


3573 


194 


8912(243) 


26644(313) 


8039(135) 


3414(65) 




2005-08-08" 


3590 


211 










330(33) 


2006-02-12" 


3778 


399 


2213(121) 


4960(135) 


1955(67) 


2283(54) 




2006-03-18" 


3812 


433 










260(26) 


2006-08-03« 


3950 


571 










340(34) 


2006-08-08« 


3955 


576 


590(66) 


729(53) 


867(47) 


1028(37) 




2007-02-20* 


4151 


772 


314(51) 


166(27) 


211(26) 


363(26) 




2007-03-18* 


4177 


798 










110(11) 


2007-07-16* 


4297 


918 










23(13) 


2007-07-24'' 


4305 


926 










28(14) 


2007-08-07'' 


4319 


940 




138(26) 




296(25) 




SN 2006bp 


2007-05-14' 


4234 


401 


355(53) 


257(33) 


464(36) 


618(37) 




2007-06-03' 


4254 


421 










119(15) 


2007-12-19^ 


4453 


620 










150(15) 


2007-12-27^' 


4461 


628 


295(49) 


154(27) 


391(34) 


593(35) 




2008-02-07^ 


4504 


671 










111(12) 


2008-05-13^ 


4600 


767 


285(48) 


160(27) 


445(35) 


585(36) 




2008-05-17^ 


4604 


771 










107(16) 


2008-06-05^ 


4622 


789 










130(13) 


2009-12-29* 


5194 


1361 


283(49) 


140(26) 








SN 2006my 


2007-02-17* 


4148 


205 


967(91) 


528(50) 


542(41) 


803(37) 




2007-07-03^' 


4285 


342 


787(84) 


345(41) 


520(37) 


708(33) 




2007-07-11* 


4292 


349 










101(25) 


2007-07-14^ 


4295 


352 










126(25) 


2007-07-23* 


4304 


361 










133(25) 


2008-02-06^' 


4502 


559 


731(82) 


312(40) 


504(37) 


742(35) 




2008-02-15^ 


4511 


568 










126(35) 


SN 2006OV 


2004-06-10" 


3166 


-897 


508(109) 


356(53) 


789(114) 


822(214) 




2007-07-03^ 


4284 


321 


606(108) 


415(53) 


752(112) 


1162(200) 




2008-02-06^ 


4502 


539 


567(110) 


363(53) 


780(118) 


815(221) 




SN 2007OC 


2008-07-12' 


4659 


251+ 


477(62) 


510(48) 


297(30) 


388(27) 




2008-12-23' 


4823 


415+ 


287(49) 


192(31) 


147(22) 


214(20) 




2009-07-30'" 


5042 


634+ 


140(36) 


22(12) 








2009-12-02'" 


5167 


759' 


174(39) 


46(16) 









Notes. "FID 69 Fazio et al.; ''FID 61009 Freedman et al.; "FID 3124 Alexander et al.; "^FID 3248 Meikle et al.; "FID 20256 
Meikle et al.; ^FID 484 Gorjian et al.; «FID 30292 Meikle et al.; ''FID 40410 Rieke et al.; "FID 30496 Fisher et al.; 
^FID 40619 Kotak et al; *FID 30945 Kenney et al.; 'FID 50534 Andrews et al.; "FID 60071 Andrews et al. 
+ From date of discovery. 



